Revealing dressed-quarks via the proton's charge distribution 
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The proton is arguably the most fundamental of Nature's readily detectable building blocks. 
It is at the heart of every nucleus and has never been observed to decay. It is nevertheless a 
composite object, defined by its valence-quark content: u + u -\- d - i.e., two up [u) quarks and 
one down (d) quark; and the manner by which they influence, inter alia, the distribution of charge 
and magnetisation within this bound-state. Much of novelty has recently been learnt about these 
distributions; and it now appears possible that the proton's momentum-space charge distribution 
possesses a zero. Experiments in the coming decade should answer critical questions posed by this 
and related advances; and we explain how such new information may assist in charting the origin 
and impact of key emergent phenomena within the strong interaction. Specifically, we show that 
the possible existence and location of a zero in the proton's electric form factor are a measure of 
nonperturbative features of the quark-quark interaction in the Standard Model, with particular 
sensitivity to the running of the dressed-quark mass. 

PACS numbers: 13.40.Gp, 14.20.Dh, 12.38.Aw, 12.38.Lg 
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Experiments during the last decade have imposed a 
new ideal. Namely, despite its simple valence-quark con- 
tent, the proton's internal structure is very complex, with 
marked differences between the distributions of charge 
and magnetisation. The challenge now is to explain the 
observations in terms of elemental nonperturbative fea- 
tures of the strong interaction. In this connection, we 
demonstrate herein that the behaviour of the proton's 
electric form factor in the 6 -10 GeV'^ range is particularly 
sensitive to the rate at which the dressed-quark mass runs 
from the nonperturbative into the perturbative domain 
of quantum chromodynamics (QCD), the strong interac- 
tion sector of the Standard Model. 

The proton's momentum-space charge and magnetisa- 
tion distributions are measured through combinations of 
the two Poincare-invariant elastic form factors that are 
required to express the proton's electromagnetic current: 
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where Q = p' —p, u{p) and u{p') are, respectively, spinors 
describing the incident, scattered proton, and Fi^2[Q^) 
are the proton's Dirac and Pauli form factors. The charge 
and magnetisation distributions [l| 

Ge{Q^)^F,{Q^)^tF2{Q^), (2) 
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where r — Q'^ /[Am'j,^], ttin is the proton's mass, and e is 
the polarisation of the virtual photon that mediates the 
interaction in Born approximation. 



The first data on the proton's form factors were made 
available by the experiments described in Ref. [2|. In 
Born approximation one may infer the individual contri- 
bution from each form factor to the cross section by using 
the technique of Rosenbluth separation [3| . Namely, one 
considers the reduced cross-section, a^, defined via: 
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It is plain from Eq. (|1]) that ctr is linearly dependent on e; 
and so a linear fit to the reduced cross-section, at fixed Q^ 
but a range of e values, provides G^{Q^) as the slope and 
tG'\j{Q'^) as the e = intercept. Owing to the relative 
factor of r, however, the signal for G\j{Q^) is enhanced 
with increasing momentum transfer, a fact which compli- 
cates an empirical determination of the proton's charge 
distribution for Q^ > IGeV^. Notwithstanding this, of 
necessity the method was employed exclusively until al- 
most the turn of the recent millennium and, on a domain 
that extends to 6 GeV^ , it produced 
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and hence a conclusion that the distributions of charge 
and magnetisation within the proton are approximately 
identical on this domain [J, la]. Significantly, this outcome 
is consistent with the, then popular, simple pictures of 
the proton's internal structure in which, e.g., quark or- 
bital angular momentum and correlations play little role. 
The situation changed dramatically when the combi- 
nation of high energy, current and polarisation at the 
Thomas Jefferson National Accelerator Facility enabled 
polarisation-transfer reactions to be measured [6|. In 
Born approximation, the scattering of longitudinally po- 
larised electrons results in a transfer of polarisation to 



the recoil proton with only two nonzero components: Pj^, 
perpendicular to the proton momentum in the scattering 
plane; and Pii, parallel to that momentum. The ratio 
P_l/P|| is proportional to G siQ^/ G m {Q"^) 0,I1- A se- 
ries of such experiments [6[ |9|-[l2| has determined that 
Ge{Q^)/Gm{Q'^) decreases almost linearly with Q^ and 
might become negative for Q^ > 8 GeV^. Such behaviour 
contrasts starkly with Eq. ([6]) ; and since the proton's 
magnetic form factor is reliably known on a spacelike 
domain that extends to Q^ « 30 GeV^ [13|, ll4| , the evo- 
lution of this ratio exposes novel features of the proton's 
charge distribution, as expressed in Ge{Q^)- 

An explanation of the discrepancy between the Rosen- 
bluth and polarisation transfer results for the ratio is 
currently judged to lie in two-photon-exchange correc- 
tions to the Born approximation, which affect the polar- 
isation transfer extraction of Ge{Q^)/Gm{Q^) far less 
than they do the ratio inferred via Rosenbluth separa- 
tion [l5|. The last decade has thus forced acceptance 
of a new paradigm; viz., the proton's internal structure 
must actually be very complex, with marked differences 
between the distributions of charge and magnetisation. 

Given that sixty years of experimental effort has thus 
far discovered only one hadronic form factor that displays 
a zero; namely, the Pauli form factor associated with the 
transition between the proton and its first radial exci- 
tation (the Roper resonance) , and that this feature was 
discovered just recently |16l4l8| . the chance that the pro- 
ton's electric form factor might become negative is fas- 
cinating. It is therefore worth elucidating the conditions 
under which that outcome is realisable before the zero is 
empirically either located or eliminated as a reasonable 
possibility. This is even more valuable if the appearance 
or absence of a zero is causally connected with a funda- 
mental nonperturbative feature of the Standard Model. 

Consider therefore a continuum computation of the 
proton's elastic form factors. This has been accomplished 
within the Dyson-Schwinger equation (DSE) framework 
jig ], which is distinguished by the feature that its ele- 
ments have a direct connection with QCD. 

To illustrate this point, we note that QCD's dressed- 
quark propagator has the form 



S{p) = l/[i-i-pA{p') + B{p')]. 
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where Z{p^) = l/A{p^) 



the 



wave-function renor- 
malisation function and M{p^) = i?(p^)/A(p^) is the 
renormalisation-point-invariant dressed-quark mass func- 
tion. In QCD with massless current-quarks, any finite- 
order perturbative computation yields M{p'^) = 0. How- 
ever, a nonperturbative solution of the DSE for the 
dressed-quark propagator (QCD's gap equation) predicts 
a nonzero mass function with a strong momentum depen- 
dence |20l . l2l| - a prediction confirmed by simulations 
of lattice- regularised QCD 22], so that it is now theo- 





FIG. 1. Poincare covariant Faddeev equation. ^ is the Fad- 
deev amplitude for a proton of total momentum P = Pq+Pd- 
The shaded rectangle demarcates the kernel of the Faddeev 
equation: single line, dressed-quark propagator; F, diquark 
correlation amplitude; and double line, diquark propagator. 



visible matter in the Universe therefore lies in the emer- 
gent strong-interaction phenomena of dynamical chiral 
symmetry breaking (DCSB) and confinement [23[. 

At the subnuclcar level, DCSB has far-reaching conse- 
quences for meson properties |19| and must be expected 
to impact just as heavily on baryons. To expose novel 
aspects of this, we first recall that the proton is a bound- 
state in quantum field theory. As such, its structure is 
described by a Faddeev amplitude, y, obtained from a 
Poincare-covariant Faddeev equation [2J| , which sums all 
possible quantum field theoretical exchanges and interac- 
tions that can take place between the three quarks that 
define its valence-quark content. With ^ in hand, the 
proton's elastic form factors may be computed once the 
associated electromagnetic current is determined. 

A dynamical prediction of Faddeev equation solutions 
obtained with the realistic interactions that describe the 
dressed-quark mass function, is the appearance of non- 
pointlike quark-|-quark (diquark) correlations within the 
proton [25|, |26|. Whether one exploits this feature in 
developing an approximation to the quark-ouark scatter- 
ing matrix within the Faddeev equation |27H29f , as illus- 
trated in Fig.[Tl or chooses instead to eschew the simpli- 
fication it offers, the outcome, when known, is the same 
[30[ . Notably, empirical evidence in support of the pres- 
ence of diquarks in the proton is accumulating |3ll435l | . 

For a proton described by the amplitude in Fig.[TJ 
the electromagnetic current is known [36[. The key el- 
ement in constructing that current is the dressed-quark- 
photon vertex. It is plain from a consideration of the 
Ward-Green- Takahashi identities 37-4(| and the struc- 



ture of the functions in Eq. ([7]) that the bare vertex (7^) 
is not a good approximation to the dressed vertex for 
Q^ < 2GcV^, where (as above) Q is the incoming pho- 
ton momentum. This has long been clear j4l| and recent 
years have produced a sophisticated understanding of the 
coupling between the photon and a dressed-fermion. Two 
model-independent results, which have emerged from the 



retically established that chiral symmetry is dynamically 
broken in QCD. The origin of the vast bulk of the mass of 



vast body of literature, are crucial herein [4^ |43|: the 
Ansatz described in Ref. [4l| is the unique form for the 
solution of the longitudinal Ward-Green- Takahashi iden- 
tity; and the transverse part of the dressed vertex ex- 
presses a dressed-quark anomalous magnetic moment dis- 
tribution, which is large at infrared momenta. Stated 
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FIG. 2. Upper panel. Dressed-quark mass function, a = 1 
specifies the reference form and increasing a diminishes the 
domain upon which DCSB is active. Lower panel. Response 
of (IpGe/Gm to increasing a; i.e., to an increasingly rapid 
transition between constituent- and parton-like behaviour of 
the dressed-quarks. Data are from Refs. [a. I9l-[l3| . 



simply, the photon to dresscd-quark coupling is markedly 
different from that of a pointlike Dirac fermion. 

The computation of the proton's elastic form fac- 
tors, using the elements detailed above, is exemplified 



in Refs. J27l - l29l |. We use that framework herein, with 
the dressed-quark mass-function illustrated in the upper 
panel of Fig. [51 the associated dressed-quark propagator, 
and the following dressed-quark-photon vertex: 



with q = k— p, t = k+p, and 41 1 



r]i''ik,p) = J2Mk,p)L^^{k,p)., 
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AB{k^,P^); and 



0(p2)]/[fc2 - p2]^ .^^ith A, B in Eq. ^. The second term 
in Eq. ([5]) expresses the momentum-dependent dressed- 
quark anomalous magnetic moment distribution, with 
<r = 0.4 being the modulating magnitude [43 1. 

In order to highlight a connection between DCSB and 
the <5^-dependence of proton form factors, one may intro- 
duce a damping factor, a, into the dressed-quark prop- 
agator used for all calculations in Refs. [27'|. [Explicitly, 



we write 63 -^ ab^ in Eq. (A. 19) of Ref. [27[, the effect 
of which is a modification in Eq. ([7|) that may be ap- 
proximated as B{p) -^ B{p){l + af{p))/{l + a'^f{p)), 
fip) = 2(p/2)V(l + b/2)^).] The value a = 1 speci- 
fies the reference form of the dressed-quark propagator, 
which was obtained in a fit to a diverse array of pionprop- 
erties [4J|. It produces a chiral-limit condensate 451447 
(gg)5 = — (0.250GeV =: Xo)^! ^^id is associated with a 
prediction of the pion's valence-quark distribution func- 
tion [48| that was recently verified empirically [49| . 

As a is increased, the rate at which the dressed-quark 
mass function drops towards its perturbative behaviour 
is accelerated so that, as evident in the upper panel of 
Fig. [21 the strength of DCSB is diminished and the influ- 
ence of explicit chiral symmetry breaking is exposed at 
smaller dressed-quark momenta. This is the qualitative 
impact of a that we exploit herein. 

At each value of a, we repeated all steps in the compu- 
tation detailed in Ref. 27[. Namely, we solved the Fad- 
deev equation to obtain the proton's mass and ampli- 
tude, and, using that material, constructed the current 
and computed the proton's elastic form factors. (The 
scalar and axial- vector diquark masses were held fixed as 
a was varied, in which case the nucleon mass, 771 jv, drops 
by < 1% as a is increased from 1.0 to 2.0. Since damping 
was deliberately implemented so that the pointwise evo- 
lution of M(p^) to its ultraviolet asymptote is accelerated 
without changing M(p^ = 0) and because the computed 
values of masses are primarily determined by the infrared 
value of mass- functions [50| , this is a reasonable assump- 
tion on the input and an understandable result for rriN-) 

The effect on Ge{Q^)/Gm{Q^), produced by suppress- 
ing DCSB, is displayed in the lower panel of Fig.[2l The 
impact is striking. For a = 1, one recovers the result in 
Ref. [28|, which exhibits a zero in Ge(<5^), and hence in 
the ratio, at Q^ w 8 GeV^. However, as a is increased, so 
that the strength of DCSB is damped, the zero is pushed 
to larger values of Q^, until it disappears completely at 
a ~ 2.0. Associating the curves in the upper and lower 
panels of the figure, one observes that apparently modest 
changes in the rate at which the mass function drops to- 
ward its ultraviolet asymptote have a dramatic effect on 
the location and existence of a zero in Ge{Q^)/Gm{Q^). 

In order to explain this remarkable behaviour, it is 
useful to recall Eqs. ([2]), ([3]). The magnetic form factor 
is a simple additive linear combination of the proton's 
Dirac and Pauli form factors. Therefore, small changes 
in Fi^2{Q^), arising from the differences displayed in the 
upper panel of Fig.[2land illustrated in Fig.[3l cannot have 
a large impact. On the other hand, the electric form 
factor is a difference, in which changes in the Pauli form 
factor are amplified with increasing Q^. 

Physically, the Pauli form factor is a gauge of the dis- 
tribution of magnetisation within the proton. Absent F2 , 
the proton's electromagnetic current would be like that of 
a Dirac fermion. The Q^ = value of the Pauli form fac- 
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FIG. 3. Q^-weighted proton Dirac (Q^Fi) and Pauli (Q^F2) 
form factors calculated with a = 1.0, 2.0. Q^Fi shows lit- 
tle sensitivity to the rate at which dressed-quarks make the 
transition between constituent- and parton-like behaviour and 
hence Fi, none at all. In contrast, Q^F2, which also appears 
in the definition of Ge, exhibits a measurable dependence. 



tor is the proton's anomalous magnetic moment; and the 
evolution with Q^ measures the distribution of anoma- 
lous magnetisation within the proton bound-state. In the 
DSE approach, the proton's magnetisation is carried by 
dressed-quarks and influenced by correlations amongst 
them. The latter are expressed via the Faddeev wave- 
function, obtained by reattaching the quark lines to the 
Faddeev amplitude. This wave function exhibits S-, P- 
and D-wave quark orbital angular momentum correla- 
tions in the proton's rest frame [51|. The resulting nu- 
cleon mass is 1.18 GeV, a value which accommodates the 
material negative pion-loop corrections [271, [S^ . 

Suppose for a moment that quarks arc described by a 
momentum-independent dressed- mass, as in Ref. |33| . In 
that counterpoint to QCD, the dressed-quarks produce 
hard Dirac and Pauli form factors, which yield a ratio 
I^pGe/Gm that possesses a zero at Q^ < 4GeV^. 

Alternatively, consider a proton comprised of dressed- 
quarks associated with the mass function in the upper 
panel of Fig. [5] This mass function is large at infrared 
momenta but approaches the current-quark mass as the 
momentum of the dresscd-quark increases. As we have 
explained, such is the behaviour in QCD: dressed-quarks 
are massive in the infrared but become parton-like in the 
ultraviolet, characterised thereupon by a mass function 
that is modulated by the current-quark mass. In this 
case, the proton's dressed-quarks possess constituent- 
quark-likc masses at small momenta. Thus, for all consid- 
ered values of a: these quarks possess a large anomalous 
magnetic moment at infrared momenta (in keeping with 
their large mass) [42]; -F'i,2(Q^) are insensitive to a on 
this domain; and hence so is the ratio ^pGe/Gm- 

On the other hand, as the momentum transfer grows, 
the structure of the integrands in the computation of 
the elastic form factors ensures that the dressed-quark 
mass functions are increasingly sampled within the do- 



main upon which the chiral condensate [451447 [ modulates 
the magnitude of Af (p^). This corresponds empirically to 
momentum transfers Q^ ^ 5 GeV^ . Plainly, as this chi- 
ral order parameter becomes smaller, a part of DCSB 
is suppressed, and the dressed-quarks become increas- 
ingly parton-like; viz., they are partially unclothed and 
come to behave as light fermion degrees of freedom on a 
larger momentum domain. Following in large part, then, 
from the fact that light- quarks must have a small anoma- 
lous magnetic moment 42] , the proton Pauli form factor 
generated dynamically therewith drops more rapidly to 
zero: the quark angular momentum correlations remain 
but the individual dressed-quark magnetic moments di- 
minish markedly. This is apparent in Fig.jS] 

Thus, as a consequence of suppressing the domain 
upon which DCSB is active, an effect expressed via a 
suppression of Xo iii the model used for this illustration, 
the zero in the ratio UpGe/Gm is pushed to larger values 
of Q^ , until it disappears from the currently accessible ex- 
perimental domain when Xo falls to roughly 80% of its 
unperturbed value. Indeed, in this case there is no zero 
in the computed result on Q^ > 0. 

An improvement of our study is possible; e.g., via the 
ah initio treatment of the DSEs detailed in Refs. 53ll54|. 
However, close inspection of results already obtained in 
that more sophisticated approach lends support to our 
conclusions; namely, as xJ/M(0) decreases, the proton's 
electric form factor approaches zero less rapidly. 

We explained that the fully-consistent treatment of 
a quark-quark interaction which yields dressed-quarks 
with a constant mass-function, produces a zero in 
fJ-pGEp{Q^)/GMp{Q^) at a small value of Q-^. At the 
other extreme, a theory in which the mass-function 
rapidly becomes partonic ~ namely, is very soft - pro- 
duces no zero at all. From a theoretical perspective, there 
are numerous possibilities in between. It follows that the 
possible existence and location of the zero in the ratio 
of proton elastic form factors [//pG_Ep(Q^)/GMp(Q^)] are 
a fairly direct measure of the nature of the quark-quark 
interaction in the Standard Model. They are a cumula- 
tive gauge of the momentum dependence of the interac- 
tion, the transition between the associated theory's non- 
perturbative and perturbative domains, and the width 
of that domain. Hence, in extending experimental mea- 
surements of this ratio, and thereby the proton's charge 
form factor, to larger momentum transfers; i.e., in reli- 
ably determining the proton's charge distribution, there 
is an extraordinary opportunity for a constructive dia- 
logue between experiment and theory. That feedback 
will enable substantial progress in contemporary efforts 
to reveal the character of the strongly interacting part of 
the Standard Model and its emergent phenomena. 
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